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Abstract This article presents a theoretical verification of 
the reinforced-matrix hypothesis derived from tensor equa¬ 
tions, <7 W = cr f + <7 m and £ w = £ f = £ m (Wood Sci Technol 
32:171-182, 1998; Wood Sci Technol 33:311-325, 1999; J 
Biomech Eng 124:432-440, 2002), using classical Mori- 
Tanaka theory on the micromechanics of fiber-reinforced 
materials (Acta Metall 21:571-574,1973; Micromechanics - 
dislcation and inclusions (in Japanese), pp 141-147, 1976). 
The Mori-Tanaka theory was applied to a small fragment of 
the cell wall undergoing changes in its physical state, such as 
those arising from sorption of moisture, maturation of wall 
components, or action of an external force, to obtain (cr A ) D = 
0-(<j F ) I + (1 - 0 )*(<j m ) d _ i . When the constitutive equation of 
each constituent material was applied to the equation (ct a ) d 
= (t)'((f)i + (1 - 0)*((7 m ) d _ i , the equations cr w = cr f + cr m and £ w 

r 

= £ = £ m were derived to lend support to the concept that two 
main phases, the reinforcing cellulose microfibril and the 
lignin-hemicellulose matrix, coexist in the same domain. 
The constitutive equations for the cell wall fragment were 
obtained without recourse to additional parameters such as 
Eshelby’s tensor S and Hill’s averaged concentration tensors 
A f and A m . In our previous articles, the coexistence of two 
main phases and a =o +o m and £ = £ = e m had been taken 
as our starting point to formulate the behavior of wood fiber 
with multilayered cell walls. The present article provides a 
rational explanation for both concepts. 
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Introduction 

Softwood fiber has two main parts: a thick secondary wall 
and thin compound middle lamella (CML). A middle layer 
(S2) forms the thickest layer of the secondary wall and is 
reinforced with fibrous polysaccharide crystals oriented 
more or less parallel to the fiber axis (see Fig. 1). The orienta¬ 
tion angle (microfibril angle; MFA) determines the mechani- 
calpropertiesofthesoftwoodfibers,includingthe longitudinal 
Young’s modulus, 13 anisotropic drying shrinkage, 45 surface 
growth stress, 6 8 and longitudinal tensile creep deformation. 9 
Consequently, the S2 layer exercises an important role in 
determining the macroscopic properties of clear wood speci¬ 
mens, with the mechanical interaction between the reinforc¬ 
ing polysaccharides and the encrusting matrix substance 
ultimately controlling the material properties of the wood. 

Any formulation of the material properties of a wood 
based on its structural hierarchy starts by considering the 
mechanical properties of the two-phase structure of the cell 
wall layer; that is, the reinforcing fibrous polysaccharide and 
the amorphous matrix. The reinforced-matrix hypothesis 
was originally proposed by Barber and Meylan 10 and pro¬ 
vides a theoretical description of the mechanical interaction 
between these two phases. This model can be expressed 
using tensor equations: ’ ’ 

cr w = cr f + cr m , and £ W = £ f = £ m . (1) 

A physical interpretation of the tensorial quantities cr w , cr f , 
a m , £ W , £ f , and £ m in Eq. 1 can be formulated as follows. The 
fibrous components of the polysaccharides are dispersed 
uniformly in each cell wall layer to form the framework 
fiber bundle. Similarly, the lignin-hemicellulose compound 
is diffused in each layer to provide the isotropic matrix 
skeleton. It is postulated that the framework bundle and the 
matrix skeleton occupy the same domain at the mesoscopic 
limit (= layer level). As a consequence, cr w can be regarded 
as the stress tensor in the cell wall fragment as the whole, 

r 

and a and cr m as the stress tensors in the polysaccharide 
framework bundle and matrix skeleton, respectively, with 
£ W , £ f , and £ m being the respective strains. 
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Fig. la-d. Structural hierarchy 
of wood (softwood), a Timber, 
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Previously we have used both Eq. 1 and these interpreta¬ 
tions as to the physical meanings of the tensorial quantities; 
that is, “the reinforced-matrix hypothesis,” to develop 
theories to explain the MFA-dependent properties of clear 
wood specimens. 6,11,12 However, neither Eq. 1 nor the 
physical interpretations of a , a, cr m , e , £, and £ m have 
been defined in any rational manner. This article resolves 
this issue using the classical theory of micromechanics 
developed by Eshelby 13,14 and Mori and Tanaka. 15 We then 
compare our formulation with that of Cave who proposed 
a constitutive relationship for the lignocellulose material of 
the wood cell wall, with the ultimate goal of validating the 
reinforced-matrix hypothesis as one of the basic theories of 
cell wall mechanics and physics. 


Mathematical derivation of the reinforced-matrix 
hypothesis 

Application of Mori and Tanaka’s theory to two-phase 
structure of the cell wall 

When material containing an ellipsoidal inhomogeneity is 
subjected to external stimulation such as heating, cooling, 
or external loading, a dimensional misfit can arise between 
the inhomogeneity and the matrix. This can induce an inho¬ 
mogeneous stress disturbance inside the current inhomoge¬ 
neity and in the matrix around it. Eshelby proposed an 
“Equivalent Inclusion Method” as a convenient way of cal- 

IT 1/1 IT 

culating this stress disturbance, a 00 . ’ ’ His method was 
intended to apply primarily to the case in which an inhomo¬ 
geneity exists at sufficient distance from the free surface of 
the material. However, in actual practice, any material 
has a free surface, and, furthermore, contains numerous 
inhomogeneities. This free surface and the neighboring 
inhomogeneities affect the stress field in both the current 
inhomogeneity and the surrounding matrix. This effect is 
called background stress. ’ Mori and Tanaka refined 
Eshelby’s method to include the case in which the stress in 

■i c -i r J 

the material is the background stress. ’ 

Consider a small fragment of the secondary wall that 
occupies a closed domain D. First, assume that domain D 
is surrounded by a free surface and incurs a dimensional 


change on moisture sorption or substance deposition (Frag¬ 
ment 1). Such changes in the physical state often generate 
some dimensional misfit between each inhomogeneity and 
the matrix around it. Second, consider the case in which 
domain D is subjected to an external loading (Fragment 2), 
and in Fragment 2, assume that no physical change of state 
occurs except that induced by the external loading. Finally, 
consider a small domain D in an actual cell wall that incurs 
a definite external stimulation in addition to the boundary 
force (Fragment 3). We can now calculate the stress distri¬ 
bution in Fragment 3 by superposing the stress fields in 
Fragments 1 and 2. In either case, an infinitely long, rod-like 
polysaccharide microfibril embedded in the matrix sub¬ 
stance can provide an ellipsoidal inhomogeneity with an 
infinitely large aspect ratio. 

In Fragment 1, the stress inside an arbitrary inhomoge¬ 
neity can be given as the sum of the background stress, a b , 
and Eshelby’s solution, crj 0 . 15,17 In contrast, the stress in the 
matrix is equal to the background stress, a \. 15,17 

In Fragment 2, we can assume that the external loading 
induces a stress distribution, cr A , inside the domain, where 
hypothetically no inhomogeneity is included. When a homo¬ 
geneous stress field is disturbed by an inhomogeneity, it can 
be modeled by Eshelby’s solution, a 2 . According to Mori 
and Tanaka’s theory, ’ the stress in each inhomogeneity is 
given as the superposition of the background stress, cr b , 
Eshelby’s solution, o 2 , and the externally induced stress, 
cr A . In contrast, the stress in the matrix is provided by super¬ 
position of cr b and a A . 

In Fragment 3, the stress distribution in D can be calcu¬ 
lated by superposing the stress fields in Fragments 1 and 2; 
that is, (T = cr b + o\ + ai + a 2 + ° A inside each inhomogene¬ 
ity, and (T = a\ + o\ + <r A in the matrix. If we denote o\ + 
o\ = cr b and crj° + o 2 = cr 00 , the stress distribution inside the 
small domain D becomes 


a = < 


<r b + <j°° + a A (= cf f ) (in I) 


cr b + cr A (= cr M ) (inD -1). 


.M 



where I is the domain of the inhomogeneities in D, and D-I 
is that of the matrix. The problem now resolves itself into 
establishing the constitutive relationship of the small frag¬ 
ment of the cell wall by analyzing the behavior of cr b , cr 00 , 
and o A in D. 
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The volume average of the internal stress distribution 
over domain D is given as cf a (see Appendix); that is, 

— f adV = — [ a A dV = 

D J D J 

D D 



where D is the volume of domain D. The above equation 
can be rewritten as 



—I jadV+ J odV 

D v I D-I ) 


\ 


+1 j (a b + o A )dV, 


-jt J(a b + o A + o°°)dV 



leading to the equation 

( aA ) D = H{<?~) I + {<? b ) I + {e A ) l ) + (l-M(<? b ) D _ l (4) 

+ <Od-.)> 


where 0 (= I/D) is the volume fraction of the inhomogene¬ 
ities in the small volume D, and (#) D , (#) I? (#) D _i are the 
integral averages of the tensorial quantities # in the respec¬ 
tive domains. We can identify the stress ( a , °°) I + (a b ) I + (<j a ) i 
with the average stress in the inhomogeneity denoted by 
(cj F ) I? with (cj b ) D _ I + (ct a ) d _i regarded as that in the matrix 
region denoted by (cf m ) d _!. Thus, Eq. 4 can be rewritten as 


Ud = *■ 


a F ) I + (l-0)-(a M 


D-I 



change in the physical state of the cell wall, such as that 
resulting from water sorption, the deposition and matura¬ 
tion of the cell wall substance, or external loading. As a 
result, the dimensions of the small volume D tend to change. 
The tensorial quantities cf m , e M , and t/ 1 are defined in D -1, 
and cr F , £ f , and of in I. 

A A 

In addition, we introduce tensor fields C and C defined 
in D to satisfy 




[C M (inD-I) 
[0 (ini) 



JO (inD-I) 
lC F (ini) 



The tensorial quantities C M and C F define the stiffness of 
the matrix substance and the crystalline polysaccharide 
microfibrils, respectively. These quantities do not depend 
on their positions within their respective domains. 

Using these tensors, we can formulate the generalized 
constitutive relationships of the matrix substance and poly¬ 
saccharide microfibrils as equations of the tensor fields 
defined in domain D as 


(j M = C M (£ — cc), a F = C F (£-&). (7) 

Taking an integral average for each equation over domain 
D, we obtain 




where 



According to Eshelby’s equivalent inclusion method, cr 00 
becomes uniform inside each equivalent inclusion. 13,14 It is 
impossible to know the distribution of cr b accurately. 

From Eq. 4 and following the assumption of Mori and 
Tanaka 15,17 that “mean field approximation” implies (cF b )j = 
(cr b )D_i = (cr b ) D (= (cr b )), we can obtain the formula 

(f7 b ) = -^-(f7“) I (=-fC7”) 


'I 

called the Mori-Tanaka theorem. We can now estimate 
the value of (o b ) from this theorem based on the value 

of CJ . 


Deriving the equations and their physical interpretations 


The constitutive equation of volume D can be derived on 
the basis of Eq. 4. In the past, we have used Eshelby’s tensor 
of inhomogeneity for this purpose. In this study, however, 
we avoid the use of Eshelby’s tensor and start with Eq. 4'. 

We first introduce stress fields & M and <r F , and strain 
fields £ and a, defined in D as 


(j M 



[(j M (inD-I) 
[0 (ini) 


(7 = 


0 (inD-I) 
[<J F (ini) 


e 


e M (inD-I) 
\e F (ini) ’ 


a 


a M (inD-I) 
a F (ini) ’ 



where £ M and £ F are the strain in the matrix and that in the 
polysaccharide microfibrils (= inhomogeneities), respec¬ 
tively. of 1 and of are the eigen-strains generated in the 
respective components. These are induced by a specific 


cr 


cr 


m 


= = D -v 

V ^ D J 

1 J e M dV, a D _ t e= 1 


V 

.M 


D-I 


D-I 


D-I 


J a M dV, 

D-I 


(e F )=jj £ F dV, (a F )=jja F dV, 
C m =(l-0)-C M , C f =fC F . 




/V /V r ^ ^ r 

As is clear from their definitions, CF m , a , C m , and C consti¬ 
tute fields of uniform tensors that are defined at every point 
in D. From Eq. 9 and Eq. 4', we can extract 


(7 m 


+ CF 




provided that (cf a ) d is regarded as the field of mean 
stress in cell wall fragment D and is identical to cf w in 
Eq. 1. 

From Eq. 9, <r m can be taken as the stress tensor of 
the matrix skeleton, and <7 r as that of the polysaccharide 
framework bundle. Similarly, C m and C r are the elastic 
constants of the matrix skeleton and the polysaccharide 
framework bundle, respectively. This implies that both the 
matrix skeleton and framework bundle occupy an identical 
volume at the mesoscopic limit. This corresponds to “the 
idea of the coexistence of two main phases in the cell 
wall.” 

Furthermore, the two formulas of Eq. 8 can be used 
to model the relationships between elastic stress and 
stiffness in the matrix skeleton and the polysaccharide 
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framework bundle, respectively, and these two formulas 
should be regarded as defining the constitutive relationships 
between the matrix skeleton and that of the polysaccharide 
framework bundle. From this perspective, we can deduce 
that 


( eM )D-I-(“ M )D-I = ^ m -“ m ’ ( eF )l-(« F )l = ( 12 ) 

where £ m is a tensor field uniformly defined at every point 
in D, and is considered as the strain of the matrix skeleton, 
and £ r is that of the polysaccharide framework bundle; a m 
and a 1 are their respective eigen-strains. Because these also 
provide the fields of the tensors uniformly defined in D, Eq. 
8 becomes 


(j m 



a 1 = C f (€*-&'). 



With the behavior of the matrix skeleton on deformation 
being entirely consistent with that of the framework bundle, 
given that the matrix skeleton and framework bundle 
coexist in identical domain D, we obtain 



Hence, the Mori-Tanaka theory provides both a theoretical 
description and a physical interpretation of Eq. I. 15 It can 
be noted that <r m , <r f , C m , C r , £ m , £ r , a m , and a L are the tensor 
fields uniformly distributed in small domain D; they were 
denoted in our previous articles by cr m , cr f , C m , C f , £ m , e\ oT, 
and a , respectively. ’ ’ ’ Consequently, we can obtain a 
- <7 f + a m from Eq. 11 and £ w = £ r = £ m from Eq. 13, which 
are the formulas of Eq. 1. 

From Eqs. 11, 8', and 13, we can obtain the constitutive 
relationship of cell wall fragment D as 

cr w = (C m + C f )£ w - (C m a m + CW). (14) 

This provides the starting point for analyzing the mechani- 

/" n -i i -| r\ 

cal behaviors of the multilayered wood fibers. ’ ’ ’ 


Simulating the behavior of the wood fiber using Eq. 14 

The strains a L and a m generated in each layer of a wood cell 
wall tend to cause dimensional changes in each wood fiber 
along its length (£ L ) and diameter (£ x ). The basic formula 

S' r-j -\ -i ^ X 

that describes such strains can be derived from Eq. 14. ” ’ 
The dimensional changes in an isolated wood fiber can be 
simulated by using this basic formula, but with the time- 
dependent patterns of a [ and a m optimized to obtain rea¬ 
sonable values of e L and £ T that are quantitatively compatible 
with observed phenomena. We can estimate the micro¬ 
scopic behavior of the constituent materials on the basis of 
optimized values of cc L and a m , and hence provide one of the 
most important purposes of the simulation. 12 

Strictly speaking, a 1 and a m are the eigen-strains of the 
polysaccharide framework bundle and that of the matrix 
skeleton, respectively, and therefore do not necessarily 
reflect intrinsic properties of the constituent materials. In 
contrast, ((/%_! and ((/\ give direct information on the 
microscopic properties of the constituent materials. The 
values of ((/%_! and ((/\ can be estimated through simula¬ 
tion by first optimizing the relationship between the macro¬ 


scopic strain, £ w (= £ m = £ f ), and the eigen-strains, cc m and 
a f . The behaviors of the microscopic strains, (£ m ) d _i and 
(£ f ) i? can be ascertained by instrumental analysis such as 
X-ray diffraction. From Eq. 12, we can determine the 
values of (a M ) D _ I and which are the microscopic eigen- 
strains caused in the polysaccharide microfibril and the 
matrix, respectively. 


Comparison with Cave’s formulation 

Outline of Cave’s formulation 

Cave is one of the pioneer researchers in cell wall micro¬ 
mechanics. He started his formulation with Eq. 4' 5,16 as 
hypothesized by Hill: 19 

(l-0)-(<7 M ) D _I + f(<7 F )l = <7 W . (15) 

Cave then introduced the constitutive relationships of two 
constituent materials in a shrinking cell wall fragment as ’ 

(<0„=c“(( £ *v,-(«“}„), 

(<A=c'« £ V(« p },)- <16) 

Equations 15 and 16 then yield 

cj w = (1-0)-C m (£ m ) d l + ^C F (s F ) l (17) 

-(1 -0)-C M (a M ) D _ I -0-C F (a F ) I . 

To develop this formula into the constitutive relationship 
of the cell wall fragment as a whole, it is necessary to relate 
the microscopic strains, (£ m ) d _i and (£ F ) I? to the macroscopic, 
£ W . Hill had introduced the following conditions on the 

1 Q 

basis of a representative volume element: 

(e M ) D _j= A M f w , {e*)=A*e\ ( 18 ) 

where A M and A F are called the averaged concentration 
tensors and depend on the elastic constants, shapes, and 
contents of the components. By using this subsidiary condi¬ 
tion, Cave developed Eq. 17 to provide the constitutive 
relationship, 

cr w = [(i - 0)-C M A M + 0 -C f A f ]e - (i - 0)-C m (a M ) D , (19) 

-fC F (a F ) r 

In his original report, Cave postulated </ = 0 by assuming 
that highly crystallized cellulose does not react with water 
molecules. 


Comparison between Cave’s formulation and the present 
one 

To adopt Eq. 19 as the constitutive relationship for simula¬ 
tion, we need to provide concrete values for A M and A F . It 
is nontrivial to determine either, because little information 
exists on the microscopic topology of the two-phase 
structure. Hence, we need to assume temporizing values 
for these factors. Cave considered that A M and A F to be 
of minor importance by assuming, along with W. Voigt 
(1889), that strain is uniformly distributed in a two-phase 
material. 5,20 By adopting Voigt’s assumption that A M = A F 
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= I (unit tensor ), 20 Cave obtained the constitutive 
relationship, 

cr w = [(1 - 0 ) • C M + 0 • C F ] £ w - [(1 - 0) • C M (a M } DI (20) 

+ 0 -C F (a F ) I ], 

and used this equation for simulating drying shrinkage and 
elastic deformation of multilayered double cell walls. 

In Eq. 20, the eigen-strains, ((/%_! and are unknown 
factors with behaviors that need to be determined from the 
simulation. However, it does not necessarily follow that the 
values of and ((/\ in Eq. 20 will reflect the intrinsic 

behaviors of the constituent materials in a wood cell wall. 
Voigt’s assumption of A M = A F = I places an unrealistic 
restriction on the situation in which the elastic moduli differ 
between the inhomogeneity and the matrix. In contrast, the 
constitutive relationship (Eq. 14) describes that between 
macroscopic strain and stress without assuming any subsid¬ 
iary conditions. By using Eq. 14 for the simulation, we can 
determine both values and behaviors of ((/%_! and (of)i. 

In our previous articles, the coexistence of the two main 
phases and Eq. 1 have been taken as basic premises with 
which to model the behavior of wood fiber . 6,7,11,12 In the 
present article, we offer a rational explanation for both 
concepts. 
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Appendix 


Volume average of the internal stress distribution over a 
small domain D as given by an externally induced stress, 
<ct a ) d 

By applying Gauss’ theorem for divergence and the 
method of integration by parts to stress field, cr, in closed 
domain, D, the volume average of a in D can be calculated 
as 


1_ 

D 


\ dij dV 


1_ 

D 

1_ 

D 


\o ± 5 ]k dV = ^-\o lk 

D U D 

j <7ik*j n k dS - J <7 ikk 


a*k 


■dV 


x s dV 


l_9D 


D 



1 , a ik n k must be nil on 3D because no external force acts on 
3D (i.e., Cauchy’s boundary condition). Fragment 2 is sub¬ 
jected to an external force through 3D and therefore o ik n k 
is nontrivial on 3D. Consequently the integral on the right 
side of Eq. 21 becomes 


1 |* 

— J n^dS, 



where (a 2 ) ik [= (cr^k + + (<7 A ) ik ] is the internal stress 

in Fragment 2. In addition, both Eshelby’s solution a 2 and 
the background stress a\ generated in finite domain D must 
satisfy Cauchy’s boundary condition for 3D: [(crj 3 )^ + (cf 6 )^ 
_ q 13,15 Therefore, (ci 2 ) ik /r k in Eq. 22 should be replaced by 
( o A ) ik n k to become 


where is Kronecker’s symbol, 3D is the boundary surface 
of the domain D, and n - (n k ) is the normal vector on 3D. 
In the case when no body force acts on D, a must satisfy 
the equilibrium condition in D: V • a [= CF ikk = (dofi/dx^ + 
(3cr i 2 /3x 2 ) + (3 073 / 3 * 3 )] = 0 , where <J ik n k is the boundary force 
that balances the external force acting on 3D. In Fragment 


ijKWV 5 - ( 23 ) 

U DD 

By applying Gauss’ theorem for divergence and the method 
of integration by parts to stress field c A , the above integral 
resolves itself as (<J A ). Hence, (a) = (c A ). 







